Elementary Electronic Excitations in Graphene Nanoribbons. 
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We analyze the collective mode spectrum of graphene nanoribbons within the random phase 
approximation. In the undoped case, only metallic armchair nanoribbons support a propagating 
plasmon mode. Landau damping of this mode is shown to be suppressed through the chirality of 
the single particle wavefunctions. We argue that undoped zigzag nanoribbons should not support 
plasmon excitations because of a broad continuum of particle-hole excitations associated with surface 
states, into which collective modes may decay. Doped nanoribbons have properties similar to those of 
semiconductor nanowires, including a plasmon mode dispersing as q%J— In qW and a static dielectric 
response that is divergent at q = 2fc_F- 

PACS numbers: 



I. INTRODUCTION 

Recent improvements in the processing of graphite 
have made possible the isolation of single graphene layers 
These carbon sheets are very promising for micro- 
electronic applications because they support the electric 
field effect: by applying a gate voltage a high mobility 
two dimensional electron or hole gas can be created 
The graphene sheets could be further processed into mul- 
titerminal devices with graphene nanoribbons acting as 
wires connecting different nanodevices. This possibility 
has motivated much recent work on carbon-based nano- 
electronics. 

The electronic properties of graphene nanoribbons de- 
pend strongly on their size and geometry 0, [Ej- This 
dependence occurs in part because edges of the ribbons 
parallel to different symmetry directions require differ- 
ent boundary conditions. For nanoribbons with zigzag 
edges ("zigzag ribbons") the boundary conditions lead 
to a particle-like and a hole-like band with evanescent 
wavefunctions confined to the surfaces [5(, which con- 
tinuously evolve into the well-known zero energy surface 
states [g, 0] as the nanoribbon width gets larger. In the 
case of armchair edge nanoribbons ( "armchair ribbons" ) 
the band structure is metallic when the width of the sam- 
ple in lattice constants units has the form 3M+1, with 
M an integer, and insulating otherwise. 

With respect to the lowest energy states, graphene 
nanoribbons are one dimensional (ID) systems. However, 
their band structure is rather different than that of more 
conventional ID systems. In armchair metallic nanorib- 
bons the dispersions of the electron and hole bands are 
linear in momentum, whereas in armchair semiconduc- 
tor nanoribbons the low energy dispersion is quadratic. 
On the other hand, zigzag nanoribbons have a peak in 
the density of states at the Fermi energy. This affects 
the low energy physics and the screening properties of 
graphene nanoribbons. In this work we compute, within 
the Random Phase Approximation (RPA) the dy- 
namical polarizability and the dielectric constant of a 



doped graphene nanoribbon. We calculate the low en- 
ergy charged collective excitations, and analyze the dif- 
ferences between conventional semiconductor-based one 
dimensional systems and graphene nanoribbons. 

We now summarize our results. Within the RPA, we 
find that for most undoped armchair nanoribbons there 
are no plasmon modes associated with interband transi- 
tions. The only exception is metallic nanoribbons, which 
support a plasmon that disperses as qy '— lng, with q 
the momentum of the excitation along the nanoribbon. 
Unlike plasmons in semiconductor based ID systems, in 
metallic nanoribbons the energy of the plasmons does not 
depend on the level of doping of the system. The case 
of undoped zigzag nanoribbons turns out to be consider- 
ably more complicated, and we do not carry through a 
full RPA analysis as in the armchair chase. However, we 
argue that a broad particle-hole continuum in the zigzag 
case implies that any poles in the inverse dielectric func- 
tion will be damped, so that there cannot be propagating 
plasmon modes. Finally, we show that doped semicon- 
ductor nanoribbons support low energy plasmons that 
disperse as q\J— In q and, in the armchair case, depend 
on the density of carriers in the ribbon in a standard way. 

This paper is organized as follows. In Section II we 
introduce the Dirac Hamiltonian used for describing the 
electronic properties of graphene. In subsections II-A 
and II-B we review the non-interacting electronic struc- 
ture of armchair and zigzag nanoribbons. Section III 
is dedicated to defining a generalized dielectric function 
for graphene nanoribbons, and to presenting results for 
collective charge excitations in intrinsic (undoped) and 
doped graphene nanoribbons which may be derived from 
it. We conclude in Section IV with a summary of the 
results. 



II. MODEL HAMILTONIAN. 

Graphene is a honeycomb structure of covalently 
bonded carbon atoms, which should be treated as a tri- 
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angular lattice with two basis atoms, denoted by A and 
B, per unit cell. In graphene the low energy dispersions 
of the electron and hole bands are linear near the points 
K = —{\,-7s) and K' = of the Brillouin 

a-o v 3' a V 3' y/3 > 

zone. Here ao is the triangular lattice parameter of the 
graphene structure. In the k p approximation the 
wavefunctions are expressed in terms of envelope func- 
tions [^(r), -0s (r)] and [^(r), "0s ( r )] f° r states near the 
K and K' points, respectively. The envelope wavefunc- 
tions may be combined into a four vector which satisfies 
the Dirac equation Hty = eVP, with 
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where 7 = \fitj2, with t the hopping amplitude between 
nearest neighbor carbon atoms. Note that k denotes the 
separation in reciprocal space of the wavevector from the 
K (K 7 ) point in the upper left (lower right) block of the 
Hamiltonian. 

The bulk solutions of Hamiltonian (fT]) have en- 
ergies £=S7do|k| with s=±l and are degenerate 
in the valley index. The wavefunctions take the 



form [e- 4ek ,s,0,0]e ikr /C 
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valley, 
Here ( 



and 



arctan k x /k v , and C = ^J2L X L V is a normalization con- 
stant, with L x and L y the sample dimensions. 



A. Armchair nanoribbon. 

In an armchair nanoribbon the termination at the 
edges consists of a line of A-B dimers, so that the wave- 
function amplitude should vanish on both sublattices at 
the extremes, x=0 and x=W + clq/2, of the nanoribbon. 
To satisfy this boundary condition we must admix valleys 
■5] , and the confined wavefunctions have the form, 
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with energies e n ^(ky)=sa^^jk\ + k y . The allowed val- 
ues of k n satisfy the condition [f| 



k n — 
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For a width of the form W=(3M + l)ao, the allowed val- 
ues of k n , k n — ^ ( 2 2M+1" ) ' create doubly degenerate 
states for n 7^ — 2M — 1, and allow a zero energy state 
when k y — > 0. Nanoribbons of widths that are not of 
this form have nondegenerate states and do not include 
a zero energy mode. Thus these nanoribbons are band 
insulators. The allowed values of k n are independent of 



the momentum k v , and the coupling between the motion 
in the x and y directions only enters in the direction of 
the isospin of the state, 9k n ,k y - 



B. Zigzag nanoribbon. 

The boundary conditions for zigzag nanoribbons are 
vanishing of the wavefunction on the A sublattice at one 
edge, x=0, and on the B sublattice at the other, x=W. 
The geometry of the zigzag nanoribbon does not mix the 
valleys, and solutions near the K valley with wavevector 
k y are degenerate with solutions near the K' valley with 
wavevector -k y . 

For the K valley and a given value of k y the nanorib- 
bons wavefunctions take the form 
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where C is the appropriate normalization constant. The 



corresponding eigenvalues are e — sjao^Jky — z\ with 

s=±l. The boundary conditions lead to a transcendental 
equation for the allowed values of z n , 
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For k y > 1/W, equation ([5]) has solutions with real val- 
ues of z which correspond to states localized at the sur- 
face. There are two branches of such surface states, one 
with positive energies, and the other with negative en- 
ergies. These energies approach zero as k y increases. In 
undoped graphene the Fermi energy is located between 
these surface state branches. For k y < 0, there are no so- 
lutions with real z and surface states are absent. The 
solutions to Eq. [5] have purely imaginary z and rep- 
resent confined states. For values of k y in the range 
< k y < 1/W, surface states from the two edges are 
strongly mixed and are indistinguishable from confined 
states. In zigzag nanoribbons the allowed values of z n 
and the corresponding wavefunctions in the transverse 
direction x depend strongly on the value of the longitu- 
dinal momentum of the state k y . This is a consequence 
of the chirality implicit in the Dirac Hamiltonian, Eq.([T]). 



III. DIELECTRIC CONSTANT. 

In nanoribbons the electrons are confined to move just 
in the y direction, and the momentum in this direction, 
k y , is a good quantum number. The motion in the x di- 
rection is confined and it is described by a band index n 
and the conduction or valence character s. The form of 
the dielectric constant depends strongly on whether the 
motion in the longitudinal and transverse directions of 
the wire are coupled, as in the case of the zigzag nanorib- 
bon, or decoupled, as in the armchair nanoribbon. In 
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the latter case the dielectric constant can be described 
as a tensor in the subband indices whereas in the for- 
mer case it also depends on the momentum k y of the 
electronic states. For undoped ribbons, we will analyze 
the armchair case in detail, whereas for the zigzag case 
we will argue that propagating plasmons are unlikely to 
exist. The results for doped graphene nanoribbons are 
qualitatively similar to those of standard semiconductor 
nanowires for both cases. 



A. Armchair nanoribbon. 

For armchair nanoribbons the motion of the carri- 
ers on the wire is coupled with the transverse motion 
only through the relative phases of the wavefunctions [cf. 
Eq.@_L I n this situation the generalized dielectric func- 
tion [Tl|, [ll] for the nanoribbon is given by 



-ijrnn 



(q, u>) = 5 im 8j„ - Vij mn (g)n m ,„ (q, to) (6) 



where the function Tl mn (q,uj) is a generalized irreducible 
polarizability function. In this expression i,j,m,n are 
composite indices specifying the subband index and 
the valence/conduction band character s = ±1. The 
wavevector q is along the direction in which the electron 
motion is free, and Si m is a Kronecker delta function. 
The explicit expression for II takes the form 



H n ,n>(q, u) 



E 



f{en'{k y + q)) - f(e n {k y )) ^ 
e n '{k y + q) - e n (k y ) - ^ 

Fn,n'(ky,k y + q) , (7) 



where g s is the spin degeneracy, /(e) is the Fermi distri- 
bution function, e n (k y ) is the subband dispersion, and 
F n ,n'(k y , ky + q) is the square of the overlap between 
wavefunctions, 



F„ 



f(ky,k y + q) = -(l + ss'coa9), (8) 



with 9 the angle between the wavevectors (k n ,k y ) and 
(k n i , k y + q) . The Coulomb matrix elements for the arm- 
chair nanoribbon take the form 

v ijm n(q) = v\i^ji\ m ^ n \(qW) = du du' 

Jo Jo 

cos{ir(i—j)u] cos[7r(m — n)u']v(qW\u — u'\) , (9) 

where the one-dimensional Fourier transform of the the 
Coulomb interaction has the form [l3| 

2e 2 

v g (x-x') = v(q\x-x'\) = K (q\x - x')\) . (10) 

£0 

Here eo is the background dielectric constant and Kq(x) 
is the zeroth-order modified Bessel function of the second 
kind, which diverges as Tn(:r) when x goes to zero. 

In order to obtain the last expression in Eq. [5] we have 
assumed that the Coulomb interaction neither produces 
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FIG. 1: (Color online) Different matrix elements of the 
Coulomb interaction as function of qW for a armchair ter- 
minated nanoribbon, Equation [9] 



intervalley scattering nor changes the sublattice index, 
as should be appropriate in the long wavelength limit. 
Note that the interaction is strictly zero for arbitrary q 
if \i ~ j \ + \ m ~ n \ is an odd integer. 

In Fig. Q] we plot different matrix elements of the 
Coulomb interaction as a function of the momentum q. 
Clearly vq^ is much larger than the others, and it diverges 
as — lii(qW) when q — > 0. Thus, in graphene nanorib- 
bons, both in the undoped and doped cases, the only 
relevant transitions at long wavelengths involve subband 
indices which satisfy the condition i — j=m — n=0. For 
the low energy excitations in the long wavelength limit 
we need only include the Coulomb matrix element wo,o, 
and can to a good approximation ignore the others. 

The charge density excitations are obtained by find- 
ing the zeros of the real part of the determinant of the 
dielectric matrix Eq.©. These excitations are well de- 
fined provided they are not damped by the continuum of 
electron- hole pair excitations, defined by a finite imagi- 
nary part of the determinant of the dielectric matrix. 

In undoped graphene nanoribbons at zero tempera- 
ture, only transitions from the valence subbands, s=-l 
to conduction subbands , s=l, contribute to the polar- 
izability. Moreover, if we include only vq,o, the lowest 
energy transitions are from the valence (s — — 1) to con- 
duction (s = 1) bands which have the same subband 
index n. The dielectric tensor has the form 

£(n 1 ,s),(n 1 ,-s),(n 2 ,s'),(n 2 ,-s')(<l, u) = 
u ss' u ni ,712 

The zeros of the determinant of this matrix are the zeros 
of the function 

1 - «0,0 ^( n «- ™+ + n "+.-"-) = 1 - «0,0 ^2 Xn, (12) 
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FIG. 2: (Color on/me)Real and imaginary part of the di- 
electric constant, Eq. (|12|l . versus ui for an undoped metallic 
armchair nanoribbon of width W=25ao- Dotted and dashed 
lines correspond to wavevectors gao=0.005 and gao=0.01 re- 
spectively. 




FIG. 3: (Color online) Schematic representation of the band 
structure and low energy electron-hole pair excitations in 
intrinsic (a) metallic armchair nanoribbons and (b)zigzag 
nanoribbons. In metallic armchair nanoribbons the chirality 
of the wavefunctions allows only transitions between carriers 
moving in the same direction to enter the polarizability; con- 
tributions from backscattering vanish. For a given wavevector 
q, because of the linear dispersion, all the electron hole pairs 
entering the polarizability have the same energy yaoq. In the 
case of zigzag nanoribbons the existence of surface states leads 
to a continuum in energies for electron-hole pairs with a given 
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F n - tn+ (ky,ky + q) (13) 



where A n (k v , q) = \e n (k y + q)\ + \e n (k y )\. 

In Fig. [J] we show the the real and imaginary parts 
of the dielectric constant as function of frequency for an 
armchair nanoribbon of width W—25ao- The lowest sub- 
band disperses linearly from zero so that the nanoribbon 
is metallic. Above this, the next lowest energy subband 
occurs for k\ k, Q.lia^ 1 . In FigfJ we see that the real 
part of the dielectric constant is close to unity and the 
imaginary part has a non- vanishing value for frequencies 
larger than twice the confinement e nergy o f the first sub- 
band. The frequency Q — 7ao(\/&i + <7 2 + ^i) defines 
the edge for the continuum of nonintcracting intersub- 
band electron-hole excitations. Note that the dielectric 
constant does not cross zero in the vicinity of these in- 
tersubband excitation energies. The absence of such ze- 
ros indicates there are no collective excitations associated 
with interband transitions, which in any case for energies 
greater than Kb would have a finite lifetime due to Lan- 
dau damping. 

For much lower frequencies, the imaginary part of the 
dielectric function remains zero down to an energy ~ 
yaoq, where a peak reveals the existence of electron-hole 
excitations at a single frequency. This corresponds to an 
intraband excitation in the metallic subband. (We note 
that in the lower panel of Fig. [2]this is slightly broadened 
due to the small imaginary iS added to uj to make the 



peak visible. In the limit £ — > this peak becomes a 
delta function.) 

The collapse of the particle-hole continuum into a sin- 
gle line in the w vs. q plane is very common in one dimen- 
sional metallic systems. However, its presence here turns 
out to depend in a fundamental way on the chirality of 
the Dirac fermions. This enters through the overlap F 
that appears in the polarizability, and guarantees that 
for a given q only particle-hole excitations of a single 
frequency contribute to Eq. [Tj|l i.e., A ra= o is indepen- 
dent of k y when F ^ 0. The overlap F = (1 + cos#)/2, 
with 8 the angle between (0,k y ) and (0, k y + q), is dif- 
ferent than zero only for 8 = n. For excitations across 
the Dirac point, this corresponds to forward scattering, 
with a single energy for a given q. Although there is a 
continuum of energies for a backscattered particle, these 
do not contribute because the wavefunction overlap F 
strictly vanishes. This is illustrated in Fig. [SUa). 

From the top panel in Fig. [5] one may see that the 
real part of the dielectric function vanishes at a fre- 
quency close to but just above the particle-hole excita- 
tion frequency. Since the imaginary part of the dielectric 
function vanishes at this frequency, one obtains a single, 
sharply defined plasmon mode. This is the only collective 
excitation that occurs in undoped graphene nanoribbons. 

In FigU]we plot the intraband plasmon frequency as a 
function of qW . It is clear in the figure that in the long- 
wavelength limit, qW < 0.001, the Bessel function can 
be approximated by -ln(qW) and the plasmon disperses 
as qy/— \n(qW). Keeping only intraband transitions, the 
effective polarizability entering in Eq. [T^] for metallic 
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FIG. 4: (Color online) The continuous line indicates the plas- 
mon dispersion in a metallic graphene nanoribbon calculated 
within the RPA. The dashed line shows the local longwave- 
length plasmon dispersion, Eq |16l The dotted line is the single 
particle excitation. 



armchair nanoribbons can be obtained analytically for 
finite temperatures and dopings, 



Xo(s,w,/3,m) 
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with (3 = l/ksT and /x the chemical potential. This 
result applies both for doped and undoped nanoribbons. 
In the longwavelength (q — > 0) limit we can easily show 
that the plasmon frequency is given by, 



ire 



\i/2 

y/-]n(qW)y/fi(q,l3,n) 



(16) 

Some comments on this result are in order. First, the 
plasmon has exactly the same dispersion as the normal 
one-dimensional plasmon of a metallic nanowire. How- 
ever, at zero temperature, in graphene nanoribbons the 
plasmon frequency is independent of the chemical poten- 
tial and therefore of the density. A classical one dimen- 
sional plasmon has a dependence ^/hib, where nn> is the 
one-dimensional carrier density. A similar anomaly oc- 
curs in two dimensional doped graphene where the plas- 
mon frequency behaves as (n^p ) , compared with the 
("2_d) 1 ^ 2 classical dependence Ha, . This anomaly 
in the density dependence of the plasmons is a direct con- 
sequence of the quantum relativistic nature of graphene 
[16j . Second, in undoped two dimensional graphene a 
finite temperature produces thermo-plasmon excitations 
associated with transitions between thermally activated 
electrons and holes 0, EH- In metallic nanoribbons 
the absence of backscattering precludes the existence of 
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FIG. 5: (Color online)(a) The static polarizability as func- 
tion of the wavevector of a doped semiconductor armchair 
nanoribbon for different densities. (b)Intrasubband plasmon 
dispersion of the same system. 



such thermo-plasmon excitations. The temperature only 
slightly alters the plasmon dispersion. Finally, we stress 
the strong dependence of the plasmon dispersion on the 
nanoribbon width. This indicates that plasmon spec- 
troscopy can be used as a characterization of the width 
of the nanoribbon. 



B. Doped semiconductor armchair nanoribbon. 

We next analyze the low energy excitations of an 
electron-doped armchair nanoribbon. We confine our- 
selves to the case that the Fermi energy lies in the low- 
est energy subband. In this situation the plasmon cor- 
responds to excitations within the same subband. The 
polarizability can be approximated by 



2 f F A(k ,q) 

n ,o = - / dk y F 0+fi+ (k v , k y + q) — — . 

7T J-k P [n^r 



A(k y ,qf 
(17) 

with A(k y , q) — \ei(k y + q)\-\ei(k y )\, and the collective ex- 
itations are obtained from the equation 



1 - "o,o (qOIIco (<?,w) = 



(18) 



The Fermi wavevector appearing in Eq |17l is related to 
the one dimensional density of extra carriers in the 
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nanoribbon by the expression kp = 7rni£)/2. In Fig. 02b) 
we plot the plasmon energy for this case and for different 
densities. As expected in one dimensional systems, the 
plasmon behaves as qy/— IngW. However, in contrast to 
the metallic case, the low energy subband dispersion is 
parabolic, and for small q the overlap between electron 
and hole wavefunctions, F 0+:0+ (k y} k y + q)=(l + cos 0)/2, 
is close to unity. Thus these plasmons have the usual 
•JniB dependence on the electron density. The possibil- 
ity of backscattering is also reflected in the static polar- 
izability. In Fig. 03a) we plot the static polarizability 
as a function of the momentum. This quantity is loga- 
rithmically divergent at q=2kp as in the standard one 
dimensional electron gas. This is a consequence of the 
perfect nesting and the possibility of backscattering from 
-kp to kp. It suggests that at sufficiently low tempera- 
ture, the system may be unstable to the formation of a 
charge density wave when coupling to lattice phonons is 
included. 



unity only for transitions from valence states, s=-l to 
the conduction states s=l. Unlike the armchair case, in 
zigzag nanoribbons there is a broad particle-hole contin- 
uum of excitations, from zero energy up to ~ 7, for small 
q. This occurs because of the existence of surface states 
which are essentially at zero energy over a broad range 
of k y , allowing arbitrarily low energy particle- hole exci- 
tations for any q smaller than an appreciable fraction of 
the Brillouin zone size, as illustrated in Fig. 05b). These 
necessarily damp any plasmon poles that might otherwise 
appear in the inverse dielectric function, and we expect 
no propagating plasmon modes for the undoped zigzag 
nanoribbon. 

The doped case is essentially the same as that of the 
doped armchair nanoribbon and standard semiconduct- 
ing nanowires: one may linearize the single particle spec- 
trum around the Fermi energy, leading to a narrow line 
of particle- hole excitations, above which a plasmon pole 
should appear with dispersion qy/— InqW. 



C. Zigzag nanoribbon. 

The calculation of plasmon poles in the zigzag nanorib- 
bon case is complicated by the fact that the parameter 
z n in the wavefunctions in Eq. OS is a function of k y . In 
this case the Coulomb matrix element cannot be writ- 
ten only as a function of the momentum transfer q. The 
generalized dielectric function within the RPA obeys 

^ijmniky > ky , q,LO^ — &im3jn&k y ,k y ^ijmn{ky, ky, (/) X 

f{tm{ky + q)) - f(e n (ky)) 
" £ m (k y + q) - e n (k y ) - K(u> + iS) 

Here g v =2 is the valley degeneracy. The matrix element 
of the Coulomb interaction has the form 

Vi jm „(k y ,k y ;q)= dx dx'^* ky (x) ■ $ jtky+q (x) x 
Jo Jo 

v q (x - x')y* m}k , y+q (x') ■ ^ nK {x'\20) 

where vj/* . vj/ denotes a dot product of the vectors ap- 
pearing in Eq. 03 The collective excitation spectrum is 
obtained by the condition that the determinant of the 
dielectric matrix, Eq. (fT9)) , vanishes. In the case of un- 
doped graphene the dielectric function is different from 



IV. SUMMARY 

In this paper we have analyzed the collective mode 
spectrum of graphene nanoribbons. In the undoped case, 
we find only metallic armchair nanoribbons support a 
propagating plasmon mode. This mode mode may prop- 
agate undamped because the chirality of the wavefunc- 
tions prevents it from decaying into particle-hole pairs. 
We argued that undoped zigzag nanoribbons will not sup- 
port plasmon excitations because this supports a broad 
continuum of particle-hole excitations with low wavevec- 
tors, and there is nothing to protect a plasmon from de- 
caying into this. Doped nanoribbons turn out to have 
properties similar to those of semiconductor nanowires, 
including a plasmon mode dispersing as qy/ — In qW with 
a coefficient that vanishes when the doping vanishes (ex- 
cept in the metallic armchair case) , and a static dielectric 
response that is divergent at q = 2k p . 
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